Entanglement between two superconducting qubits via interaction with non-classical 

radiation 
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We propose a scheme to physically interface superconducting nano-circuits and quantum optics. 
We address the transfer of quantum information between systems having different physical natures 
and defined in Hilbert spaces of different dimensions. In particular, we investigate the transfer 
of the entanglement initially in a non-classical state of a continuous-variable system to a pair of 
superconducting charge qubits. This set-up is able to drive an initially separable state of the qubits 
into an almost pure, highly entangled state suitable for quantum information processing. 
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Control of the dynamics of a complex quantum sys- 
tem requires a trade-off between tunability and protec- 
tion against noise. To this end one can be interested in 
processes where some physical properties of a subsystem 
are reliably transferred onto the state of a second one 
(of perhaps different nature) where information can be 
manipulated. The connection between the two subsys- 
tems is effectively realized via a physical interface. An 
interface is a communication channel used to connect the 
elements of a quantum register to perform quantum in- 
formation processing or a physical mechanism that gives 
full access to the system under investigation and allows 
to manipulate it. 

To investigate this problem, in this paper we describe 
the coupling between a nano-electronic circuit imple- 
menting a pair of quantum bits and a two-mode electro- 
magnetic field. We discuss a mechanism for the transfer 
of entanglement from a two-mode squeezed state to the 
pair of qubits. Here, the information sheltered in the 
electromagnetic medium may be manipulated, using just 
single-qubit operations, when transferred to the solid- 
state subsystem. This may offer advantages with respect 
to integrability and scalability. In particular we con- 
sider a pair of (initially independent) Superconducting- 
Quantum-Interference-Devices (SQUIDs) implementing 
two charge qubits , whereas each field mode is modeled 
as an LC circuit (Fig. [TJ. The SQUIDs can be individ- 
ually addressed by gate voltages V g whereas an external 
magnetic flux <Ag~ t allows to change the Josephson cou- 
pling Ej((fi ex t) [lland to modulate the interaction among 
the subsystems 0, 0- Direct experimental evidence of 
the use of these systems as controllable coherent two- 
level systems has already been provided 0, 0| . 

We first analize a single SQUID plus LC oscillator 
(Fig. Gib)). We introduce the phase drop across the 
SQUID (tp a ) and across the LC circuit (<fib)- The con- 
jugated variables are the excess charge on the SQUID is- 
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FIG. 1: (a) Set-up for an entanglement transfer process via 
the interface between quantum correlated field modes and a 
pair of charge qubits. Each SQUID is threaded by an external 
magnetic flux to modulate Ej(<f> ex t). (b) Equivalent circuit 
for the single SQUID capacitively coupled to a single field 
mode, modeled as a LC oscillator. 



land (Q a ) and the charge on the oscillator's capacitance 
(Pi). The Hamiltonian describing the system is: 
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where C = V/(C + C c ) 
C 2 = V/C c , V = {C d 
Ej{<p ext ) = 2£°cos(2e<A 



(1) 

Ci = V/(C g + 2C Jo + C c ), 
C c )(C g + 2C Jo ) + C C C and 
ext /h). The SQUID Hamilto- 



nian H squ id can be tuned by V g and (f> e xt- The field mode, 
described by the oscillator H em , has effective frequency 
lu = (L0C2) which comes from the inductance L a and 
the total capacitance C2 seen by the charge Ph. The cou- 
pling Hamiltonian H c describes the Coulomb interaction 
between the charges Q a and P&. 

We assume large charging energy, e 2 /2C ^> Ej(4> ext ), 
and low temperatures T <C e 2 /2C. In this regime the 
SQUID can be described by the states \m) s (m = 0, 1) 
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representing m Cooper pairs in excess in the island, 
and implements a charge Josephson qubit 0. Typi- 
cal values of C., Q ~ 1CT 15 F and C g ~ 10" 17 F guar- 
antee e 2 /2C ~ liT»£$ ~ 100 miC. In this system, 
the main sources of decoherence are noise of electrostatic 
origin, voltage fluctuations of the circuit Q, or stray 
polarization due to charged impurities located close to 
the device J\. If we set C g V g = e, |0) s and |l) s have 
the same electrostatic energy and the SQUID is not af- 
fected, at first order, by this charge noise 0>0|- At this 
working point, H squid = \Ej((t> ext )a z . s with a computa- 
tional basis {|+> , where |±) s - (1/V§) (|1> ± |0» a 
are eigenstates of <r z . s splitted by Ej/h ~ 10GHz. 
We introduce the operators d and d^ (la, aJ~\ = 1) via 
P b = (huC 2 /2) 1/2 (a + at), ^ 6 = i(h/2uC 2 ) 1/2 (a - fit) 
and we get i? em = ftw (d^d + 1/2). Taking C2 — 1 pF 
and F Q ~ 10 nff , achievable by present day technology, 
we have ui ~ 10 GHz. The coupling between the SQUID 
and the field mode can be tuned on and off resonance 
by modulating the energy splitting of the qubit via 4> ex t ■ 
If Ej(<j) ext )/h is set to be much different from u>, the 
coupling is effectively turned off and the qubit evolves 
independently from the field mode. On the other hand, 
for the quasi-resonant condition Ej((j) ext ) ~ &j, we use 



Q a = lea 



x> „, with & x a 
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-) s (+|), so that 



Ml [(der + . s + dt<7_. s ) + (a) & + , s + a<T_ )S )l , (2) 



where ft = e\j2uC2/hC 2 is the Rabi frequency of the 
interaction and a+ tS = a_ s = |+).(— |- The Hamil- 
tonian Eq. @ is frequently found in quantum optics 
problems. The first and second term preserve the to- 
tal number of excitations in the system and allow for the 
restriction of the computational basis in Hsquid ® H em 
to {|— ,n) , |+, n — 1)} 8 . The other (counter-rotating) 
terms induce leakage from this subspace. They can be 
neglected in the Rotating Wave Approximation (RWA) 
that is applicable when Q <C u, Ej(cj> ex t)/fi, achieved if 
we take C c ~ 10 _17 F (weakly coupled subsystems) so 
that Ci ~ 10~ n F and ~ 0.1 GHz. In this regime the 
eigenstates of SQUID+ficld mode are entangled states 
forming a series of doublets splitted by hVly/n. It is worth 
stressing that at the working point V g = e/C g , intra- 
doublet transitions are forbidden |2|. The system is thus 
protected, to a certain extent, against decoherence. 

We now describe the interaction of a pair of SQUIDs 
with non-classical radiation. We consider the two-mode 
squeezed state \S(r)) ab = J2n=o Vn K n ) a b> where r is 
the squeezing parameter and r] n = (tanhr)™/ coshr 
The entanglement between modes a, b is a function of 
r. Squeezed microwaves can be generated off-line with 
Josephson parametric oscillators |9( and then used for our 
protocol. The SQUIDs can be integrated in the waveg- 
uides used for the transmission of the signal Js], with 
the gate-plates orthogonal to the direction of propaga- 
tion of the fields. Quality factors ~ 10 4 for a supercon- 



ducting transmission-line are within the state of the art. 
For u> ~ 10 GHz, this gives photon-lifetimes^ 1 psec , 
allowing for a coherent dynamics. The SQUIDs are pre- 
pared in a pure separable state piz(0) = Pi(0) ® pz(0)- 
The interaction between each SQUID and a field mode 
is driven by the co-rotating part of Eq. J2J, H^P wa , with 
i = 1,2 and j = a,b. The time evolution operator is 
U(t) = ®ije~i H ^ at 0|- However, the reduced state 
of the SQUIDs pi 2 (i) is inseparable, in general, be- 
cause the evolution could have transferred quantum cor- 
relations from the fields to the qubits. To see it, we 
derive the operator-sum representation of the SQUIDs 

evolution p 12 (i) = Tr ab {u(t)p 12 (0) <g> p ab (0)ft (t)\ = 

E M E™, P =o^rMo)^;r pt < where we the Kraus °p- 

erators K™ p — E^Lo^M (^n,p\U(t)\n,n^ have been 

introduced. Calculating the matrix elements of U(t) 
over the number states of the field modes, a set of five 
Kraus operators is found. If the initial state of the 
two SQUIDs is specified, a simplification is possible and 
the number of Kraus operators is reduced. We assume 
Pi2(0) = I — , — ) 12 (— , — I , that can be prepared using 
standard techniques [l|. We get the effective represen- 



tation Pl2(t) 
where 



fj,= l Z^m=0 ""-/J 



K™ = Vm cos 2 (nV^t) |— ,— >ia<— ,— I 

- r) m+ i sm 2 (nVm+ It) |+, +) 12 (-, -| , 
K™ = Vm cos(fl\/mt) sm(£l\/mt) | — , +) 12 (— > — I 
K™ = r\ m cos(fty/mt) sin(fl\fmt) |+, — ) 12 (— , — | 



(3) 



P\2(r, t) 



(4) 



K™ is responsible for zero and two-photon processes that 
leave the two field modes with the same number of pho- 
tons. K™ and K™ describe single-photon processes in 
which one of the SQUIDs absorbs an incoming photon. 
Using Eqs. 10), the density matrix of the SQUIDs, in the 
ordered basis {|ij)} 12 (i, j = +,— ), takes the form 

( A(r,t) -D(r,t)\ 

B(r,t) 

B(r,t) 

\-D(r,t) C(r,t) J 

Here A(r,t) = J2n,o Xnn(r) cos 4 (Qy/nt) , B(r, t) = 
E^o Xnn(r) sin 2 (nynt j cos 2 (Sly/fit), D(r, t) 
E^o Xnn+i(r) sin 2 (tty/nTTt) cos 2 (0^/nt) and 
C(r,t) = 1 - 2B(r,t) - A(r,t), with Xnm = i) n r\ m . 
To quantify the entanglement between the qubits, we 
choose the negativity of partial transposition (NPT). 
NPT is a necessary and sufficient condition for entangle- 
ment of any bipartite qubit state 01 ■ The corresponding 
entanglement measure is defined as Snpt = — 2A~(r, t), 
where A~ (r, t) is the unique negative eigenvalue of the 
two-qubit partially transposed density matrix /?f 2 T [llj . 
In our case, just A~(r, t) = B(r, t) — D(r, t) can be nega- 
tive for some value of r and t and it is used to compute 
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the entanglement. Enpt is shown in Fig. [2] as a function 
of the degree of squeezing r and the rescaled interaction 
time r = fit. It turns out that Enpt never becomes 




FIG. 2: Enpt versus r = fit and r. Iff Enpt > 0, there 
is entanglement between the superconducting qubits. A local 
maximum Enpt — 0-87 is achieved for r = 0.86 and f ~ 3tv/2. 

negative and only asymptotically goes to zero as r is 
increased. Once the interaction starts, the entanglement 
is transferred to the qubits, collapsing and reviving as 
the interaction time increases. The maximum of the 
transferred entanglement is £™ C P X T = 0.87, obtained for 
f ~ 37r/2 and f = 0.86 [l^. Enpt is not a monotone 
function of r as can be seen in Fig. (a) . It is known 
that the correlations in \S(r)) ab approach those of the 
maximally entangled EPR state when r — ► oo [l3j |. 
Increasing r, the contribution by higher photon-number 
terms in the squeezed state becomes more relevant. 
Each qubit is exposed to a distribution of different Rabi 
frequencies (each equal to fl^/n). These induce Rabi 
floppings at different times that interfere spoiling the 
degree of entanglement between the SQUIDs. This 
shows that a perfectly correlated continuous variable 
state can not be mapped onto a maximally entangled 
state of two qubits. On the other hand, the discreteness 
of this distribution induces the entanglement to collapse 
and revive as time goes by. This analysis is confirme d by 
considering the entanglement of formation (EoF) |14j . 
In Fig. 13 the two entanglement measures are compared, 
as functions of both r and r. From the behavior of EoF, 
we argue that almost one EPR singlet is required to 
prepare pii{f, f). 

For both the measures, the SQUIDs are separable 
just for a short initial amount of time (t ~ 1 nsec) af- 
ter which entanglement is set, persisting in time, even 
if fluctuating. The maximum of transferred entangle- 
ment is reached for t ~ 50 nsec [l5| . This is within 
the coherence time for the capacitive coupling consid- 
ered here Q and within the lifetime of the radiation, 
as already stated. A further analysis of pi2(f,f) shows 
that |B(f,f)| < \A(r,f)\,\C(r,f)\,\D(r,f)\. If, in zero- 
order approximation, we neglect B(f,t) in pi2(r,t), we 



get a density matrix close to that of the pure, non maxi- 
mally entangled state {y/A(r, t) | ) - y/C(r, t) |++))i2- 

In general, D(r,t) ^ \/ A(r, t)C(r,t), so that the state 
is mixed. The degree of mixedness in this purified ver- 
sion of pi2 (r, t) is quantified using the linearized entropy 
Si(p s (r,t)) = 4/3 [1 - Tr (pf 2 (r,t))], that ranges from 
(pure states) to 1 (maximally mixed ones). We get 
Si(pi2(r,r)) ~ 0.01 that shows that, for these parame- 
ters, the two SQUIDs are in a nearly pure state. This 
result is interesting: it has been proved, for example, 
that a bipartite mixed state becomes useless for quantum 
teleportation whenever its linearized entropy exceeds 
l-(2/[JV(JV+l)]) [3, with N the dimension of each sub- 
system. For qubits, the threshold is 2/3 ^> <S)(pi2(f, r )) 
and the state of the entangled SQUIDs could be used 
as a quantum channel in protocols for distributed quan- 
tum computation. We have calculated the purity of the 
state when B(f, f) is included, finding the same order of 
magnitude of the previous result. 

We now consider the average value of the transferred 
entanglement as the preparation of the initial state of 
the SQUIDs is varied. This allows to investigate about 
the dynamics of the superconducting qubits once dif- 
ferent separable states as (cosa |— ) + e lv sin a ® 
(cos (3 |— ) + e 1 ^ sin/3 |+)) 2 are considered. Here a, (3 € 
[0, 27r] and ip, ip £ [0, w], so that the entire Bloch sphere 
is explored. For simplicity, we take (p — ij) = and we use 
Snpt to calculate the average value of the entanglement. 
The evolution of the SQUIDs involves the complete set of 
Kraus operators. However, the density matrix pi2(r,t), 
averaged over an uniform distribution for a, /?, still keeps 
the form in Eq. Q but with more complicated matrix el- 
ements. The results are shown in Fig.^Ja). The amount 
of transferred entanglement is reduced and the peak at 
t = f , r = f is shrunk to ~ 0.4. This can be under- 
stood considering the behavior of Enpt for |+,+) 12 as 
initial state. In this case, Enpt remains negative for a 
wide range of values of r and r and has a small posi- 
tive bump for r ~ 0.6 and r ~ 1.7 (that corresponds to 

(a) (b) 

Entanglement Entanglement 
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FIG. 3: Comparison between EoF (dashed line) and Enpt 
(solid line). In (a) we plot the behavior of the two entangle- 
ment measures against the squeezing parameter r. We have 
taken r = 37r/2. In (b), the entanglement functions are plot- 
ted versus r, for r = 0.86. 
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the first peak in Fig. [2J. The previous result suggests 
that |— , — )i 2 plays a privileged role in the process of en- 
tanglement transfer. To support this idea, we look for 
the optimal preparation of the SQUIDs. Assuming small 
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entanglement-transfer from the field to the qubits can 
be efficiently tailored. This work has been supported in 
part by the UK Engineering and Physical Science Re- 
search Council grant GR/S14023/01. MP acknowledge 
IRCEP for financial support. 




FIG. 4: (a): Transferred entanglement averaged over the pos- 
sible preparations of the SQUIDs. The maximum entangle- 
ment is reduced with respect to the case of |— , — } 12 . This is 
due to the contribution from | + , +) 12 , that is separable for 
wide ranges of r and r and spoils the average entanglement, 
(b): The amount of entanglement between the qubits as a 
function of the preparation of the initial states when r — r 



and t — f. The state 



obtained for a = j3 = 0, cor- 



responds to the maximum of the transferred entanglement. 



values of a, (3, up to the their second power, we have 



pi 2 (r,t) ~N 



( a (3b + ad [id + ab — c \ 
0b + ad a' af - (3d 

ab + (3d a' f3f-ad 

\ -c af - 13d f3f~ad a" j 

with Af = (1 — a 2 — (3 2 ) and a, a! , a", b, c, d, f combi- 
nations of trigonometric functions involving both r and 
t. An analytical expression for the eigenvalues of the par- 
tial transpose is, this time, hard to be obtained. More 
insight is given by specifying r and r. In Fig.^Jb) we plot 
£npt versus a and [3 for r = r and r = f . The trans- 
ferred entanglement has a maximum equal to 0.87 for 
a — f3 — and slowly decays. This could be important, 
experimentally, because small errors in the preparation of 
the initial state do not dramatically spoil the amount of 
entanglement transferred to the qubits. The same quali- 
tative behavior is found for other values of r and t. Thus, 
the initial preparation | — ,— ) 12 provides the maximum 
achievable entanglement transfer. The entanglement of 
the SQUIDs, after the interaction, can be revealed by de- 
tecting the population B(r,r) and the coherence D(r,r) 
of the density matrix using local resonant pulses on the 
SQUIDs, along the same lines depicted in |b7f. 

We have proposed a physical interface between quan- 
tum optics and a system of two charge qubits. When 
a quantum-correlated state of light is considered, an 
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